Nonlinear Dynamics of Magnetic Islands Imbedded in Small Scale Turbulence 
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The nonlinear dynamics of magnetic tearing islands imbedded in a pressure gradient driven tur- 
bulence is investigated numerically in a reduced magnetohydrodynamic model. The study reveals 
regimes where the linear and nonlinear phases of the tearing instability are controlled by the proper- 
ties of the pressure gradient. In these regimes, the interplay between the pressure and the magnetic 
flux determines the dynamics of the saturated state. A secondary instability can occur and strongly 
modify the magnetic island dynamics by triggering a poloidal rotation. It is shown that the complex 
nonlinear interaction between the islands and turbulence is nonlocal and involves small scales. 
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Magnetic reconnection is a complex phenomenon in- 
volving plasma flows and a rearrangement of the mag- 
netic field lines inside a narrow region (the reconnection 
layer) where topologically different magnetic flux tubes 
can get interconnected and reconfigure themselves. It 
plays an important role in fusion experiments and in 
many astrophysical events [l|. In a complex fusion de- 
vice, such as a tokamak, the plasma is susceptible to 
many kinds of instabilities which can occur concurrently 
at various space and time scales. Such a coexistence of 
microturbulence and magnetohydrodynamic (MHD) ac- 
tivities has been observed in many experiments with 
some evidence of correlated effects arising from their si- 
multaneous existence. An important question to address 
is therefore the nature and amount of mutual interaction 
between microturbulence and large-scale MHD instabil- 
ities - an issue that is at the heart of multi-scale phe- 
nomena of complex systems in astrophysics, geophysics, 
nonlinear dynamics and fluid turbulence. Early ana- 
lytic attempts at investigation of this important ques- 
tion have relied on ad-hoc modeling of turbulence effects 
through anomalous transport coefficients More re- 
cently a minimal self-consistent model based on wave ki- 
netics and adiabatic theory has been used in 01 to study 
the interaction of a tearing mode with drift wave tur- 
bulence. Numerical simulation studies in Q have di- 
rectly addressed the problem of multiscale interactions 
and have taken into account the nonlinear modifications 
of the equilibrium profiles due to turbulence. Such stud- 
ies have been extended in Q to investigate the interac- 
tion between double tearing modes and micro-turbulence 
through the excitation of zonal flows. Finally in a 
numerical investigation of the interaction of a 2D elec- 
trostatic turbulence with an island whose dynamics is 
not fully self-consistent but is governed by a generalized 
Rutherford equation has been carried out. In this paper 
we report on self-consistent simulations of the multiscale 
interaction between microturbulence driven by pressure 
gradients and magnetic islands with a focus on regimes 
where the growth of the latter is essentially due to pres- 
sure effects and where small-scale dynamics appear to 
be important. The background microturbulence is found 



to induce a nonlinear rotation of the island as well as 
to signiflcantly alter its final quasi-equilibrium state by 
the excitation of a secondary instability. We discuss the 
characteristics of the various stages of the nonlinear evo- 
lution and also delineate the role of small scales in the 
overall dynamics of the system. 

We consider a minimalist two-dimensional plasma 
model based on the two fluid Braginskii equations in the 
drift approximation Q with cold ions and isothermal 
electrons. The model includes magnetic curvature effects 
and electron diamagnetic effects but neglects electron in- 
ertia and Hall effect contributions. The evolution equa- 
tions are 
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where the dynamical field quantities are the electrostatic 
potential (j), the electron pressure p and the magnetic fiux 
ip. The equilibrium quantities are a constant pressure 
gradient and a magnetic field corresponding to a Harris 
current sheet model Further, - — ~ 
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are the curvature terms with i?o representing 
the major radius of a toroidal plasma configuration. Lp 
is the gradient scale length, ta is the Alfven time based 
on a reference perpendicular length scale L± and 17^ is 
the ion cyclotron frequency. Equations ([TJS]) are nor- 
malized using the characteristic Alfven speed va and the 
length scale L^. fi is the viscosity, x± the perpendicular 
diffusivity, 77 is the plasma resistivity, = f3e/QiTA is 
the normalized electron diamagnetic drift velocity with 
/3e being the ratio between the electronic kinetic pressure 
and the magnetic pressure, p = ^ is the normalized ion 
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Figure 1: Numerical and theoretical results of the linear 
growth rate 7 versus r\ sA, v = x± = 0, A' = 6 and p — 10~^. 

sound Larmor radius. In the limit i?o 00, we recover 
the drift tearing model and when magnetic fluctua- 
tions are weak. -0^0 (ki ^ 0), the system describes 
the electrostatic interchange instability. Conversely, the 
large island limit, p = w^, = with ^iTA 1, gives the 
high (3 model which was originally introduced by H. R. 
Strauss [lo| . The minimalist model used here is, in fact, 
a reduced version of the four fields model of [11| where 
we have ignored the parallel ion dynamics and thereby 
neglected its effect on the transversal pressure balance 

As a preliminary to the numerical study of eqs.([l]|31), 
we first look at some linear results of a simplified set of 
equations with ^ = x± ^ Ki — where it is possible to 
obtain analytic relations for the linear growth rate of the 
tearing mode under the constant ip approximation 
namely, 

^ - fcjtt J_oo ^ az, 

z = -^;^x (z) + (1 + /5^a~^z2)x" (z) 

where (/> (x) = —a~^ijj (0) x (z): z — ax, a = {ilj/kyY^^ 
and A' is the standard stability parameter. The per- 
turbed pressure is given by p (x) = p'^(f>" {x) (with /5 = 
corresponding to the classical tearing situation). Fig. [1] 
shows the dependence of the linear growth rate 7(77) of 
the instability on the resistivity with p = 10~^. The 
numerical results (circles) are seen to agree quite well 
with the values (diamonds) of the solution to the an- 
alytic relation Eq.([5]). We observe that 7(77) exhibits 
a change of slope after a certain value of 77. The two 
regimes correspond to the two limiting cases Tg = or 
p — (solid line) and = (dashdot line). The in- 
tersection of the two lines gives the critical value of the 
resistivity r]c = 0.58 A'~i/2p5/2 _ g x IQ-^ for A' = 6. 
When rj > ric the linear growth rate given by the classical 
tearing case is higher than the other limiting case. The 
system chooses the more unstable case and the classical 
tearing mode is recovered with the scaling laws 7 ~ ry"^/^ 
and 6 ~ 77^/^. When rj < rjc, the coupling between p 
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Figure 2: Time evolution of the magnetic, pressure and kinetic 
energies. 



and rjj is strong and the island formation is driven by the 
pressure perturbation. The resistive layer becomes thin- 
ner or more singular than in the case Tg = where only 
a (0, t/i)— coupling exists. Further the disagreement ob- 
served between the numerical and the theoretical results 
for 77 > 5 X 10 ~^ is a consequence of the breaking down 
of the constant ip approximation in this regime. 

We now discuss the full nonlinear numerical simulation 
of Eqs.([T][3]) that explores the mutual interaction between 
small scale interchange modes and a small magnetic is- 
land. A semi-spectral code with a 2/3-dealiasing rule in 
the poloidal direction, a resolution of 128 grid points in 
the radial direction, 96 poloidal modes and that main- 
tains conservation properties of the nonlinear terms to 
a high degree, has been used. The computational box 
size is Lx = Ly = 2tt. In order to isolate the nonlinear 
mechanisms responsible for the island rotation, the lin- 
ear diamagnetic effect has been turned off in Eq. ^ . The 
effect of the latter on the evolution of the tearing mode is 
well known, namely that it leads to a real frequency and 
consequently a rotation of the island in the diamagnetic 
drift direction. Note that we have checked a posteriori 
by turning on the linear diamagnetic term in eq.(3) that 
the amount of induced nonlinear rotation (obtained by 
subtracting the linear diamagnetic frequency from the 
total rotation) remains the same. In eq.(3) we also set 
K2 = 0, since we find from our simulations that the K2 
contribution is rather weak, p and Vi, are taken to be 
equal to 1 and /?e ~ 10"^. The parameter related to the 
interchange instability is ki = 10~^. The shape of the 
equilibrium magnetic field is chosen to allow a tearing in- 
stability to develop with a poloidal mode number ky = 1 
with A' = eQ. Fig.Hshows, for p. ^ x± = V = 10"", the 
time evolution of the magnetic (E^), pressure (Ep) and 
kinetic (E^) energies of the fluctuations for r] < rjc, corre- 
sponding to a regime where the magnetic island genera- 
tion is pressure driven. Four phases are observed. First, a 
exponential growth of the magnetic island {t < ISOOr^), 
followed by a quasi-plateau phase with however an in- 
crease of the energies of the three fields (t < 4500r^). 
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Figure 4: Spectral energy densities as functions of the poloidal 
mode number ky, just (a) before {a.tt = ASOOta) and (b) after 
(at t — 5200Tyi), the bifurcation. 



Figure 3: Snapsliots of the magnetic flux and the pressure at 
t — 3000r,4 (Upper panel) and t — GOOOta (Lower panel). 



Next, a phase characterized by an abrupt growth of the 
kinetic and pressure energies in which the kinetic energy 
level equals the energy of pressure perturbations and fi- 
nally, the system reaches a new quasi-plateau phase for 
{t > 510Qta)- During the linear and first plateau phases, 
the energy associated with the pressure perturbations is 
higher than the kinetic energy, i.e the dynamics is con- 
trolled by an interplay between the magnetic flux and the 
pressure. In the second phase, t < 4500ta, the magnetic 
island is maintained by adjacent pressure cells similar to 
what is usually observed for flow cells in the nonlinear 
regime of a tearing island [J. This is illustrated in Fig. [3] 
(upper panel, t = SOOOr^). During this phase, the kinetic 
energy piles up in the flow cells which are located in the 
vicinity of the island. After t ^ 3600ta, the flow cells are 
no longer located in the vicinity of the magnetic island. 
Att> 4500r^, a sharp growth of the kinetic and pressure 
energies occurs. Far from the island the current is not 
significant, and for t/rA G [4500,5000], a dominant in- 
terchange mode outside the sheet (0ii,pii) is enhanced. 
The associated kinetic and pressure energies of the lat- 
ter are equal. Here, 4>ii means (f>{kx = l,ky ^ 1). The 
competition between the interchange and tearing modes 
lead to the generation of small scale pressure structures 
in the vicinity of the island that suffer further desta- 
bilization leading to a drastic modification of the dy- 
namics. Indeed, in less than 200 Alfven times, around 
t ^ 5000r^, an abrupt growth of the energy contained 
in the pressure perturbation is observed and the system 
dynamics chang bifurcation occurs. At larger 

times, t > 5100r^, the pressure dominates over the flow, 
Ep ^ E^, and the size of the magnetic island finally sat- 
urates. Fig. m shows the energy spectra of the fields just 
before and after the bifurcation . Before the bifurcation, 
the interchange mode is observed at fcy = 1, and as long 
as fcj^ > 2, the pressure energy is much higher than the 
kinetic energy. After this dynamical bifurcation, we ob- 
serve a persistence of small scales as well as an enhance- 
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Figure 5: Plots of the poloidal diamagnetic velocity Vdia at 
t = SOOOta (a) and t = GOOOta (b). 



ment of the energies (Fig. |3}d). We also find that a mean 
poloidal pressure and flow have been generated, and as 
we will see below, it is linked to the rotation properties 
of the island. For 8 < fcy < 50, there is a trend towards 
an equipartition of the magnetic and pressure spectra. It 
is worth noting that even though the magnetic island is 
still, at this point, in a quasilinear stage (the magnetic 
energy being concentrated on the mode ky = 1). the pres- 
sure perturbation has a fully nonlinear structure and is 
made up mainly of the modes ky < 7. An interesting 
feature clearly observed in the snapshot shown in Fig. [3] 
at t = 6000r^ is the generation of an island structure in 
the pressure field containing almost 90% of the pressure 
energy Ep. 

In the final stage where the energies reach a new quasi- 
plateau, as observed in Fig. [5J the change of dynamics 
is characterized by two important macroscopic features. 
First, there is a change of symmetry - the poloidal dia- 
magnetic velocity v^^^ ~ ^ P ha-ving even par- 
ity for t < 5100ryii (brackets mean an average over the 
poloidal direction), loses this property after the bifurca- 
tion and has in fact an odd parity in the vicinity of the 
current sheet. This change of parity is clearly shown in 
Fig. [5]) where Vdia is plotted, before and after the transi- 
tion, aXt ^ SOOOta and t = 6000r^ respectively. The sec- 
ond macroscopic change is the inversion of the poloidal 
rotation direction of the magnetic island together with 
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Figure 6: Time evolution of the poloidal position of the center 
of the magnetic island yisiandit) (solid blue line), and the 
models ydiaif) (dashed red line), yp^^{t) (dashdot green line). 

an amplification of the velocity. The amplified velocity 
arising from the nonlinear interactions is of the order of 
the linear diamagnetic velocity as verified from d posteri- 
ori runs made with the linear diamagnetic term retained 
in the equations. The change of direction in the island 
rotation can be observed in the zoomed frame of Fig. [5] 
which shows the time evolution of the poloidal position of 
the island. We find that the increase of Vdia at the transi- 
tion is linked to the coincident growth of the interchange 
mode ((/)ii,pii) which feeds the angular momentum. The 
detailed mechanism of this nonlinear geireration of an- 
gular momentum is however not known at this time aird 
remains an open question. 

Some insights into the origin of the island poloidal ro- 
tation can be obtained from Eq.([2|) where one notes that 
both the self generated zonal and diamagnetic flow terms, 
Vzon — 4> >y a^nd Vdia, Can produce a poloidal rota- 

tion of the island. To investigate the role of these flows, 
we have plotted in Fig. [5] the poloidal position of the 
center of the island y island, and the poloidal positions 
related to the contributions of the diamagnetic velocity 
Vdia and the zonal flow velocity Vzon- More precisely, we 
have plotted ydia{t) = (1/(5) J^dt l^J'^j^dx Vdia{x,t) and 

2/p,0 = (1/(5) Jgdi J%^/2^^ (""dia - Vzon)- We observe 
that the model ydia reproduces well the time evolution of 
the rotation of the island, before and after the bifurca- 
tion. At larger times, t > GOOOta , we observe that the 
contribution of the zonal flow cannot be neglected, even 
if the rotation of the island is mainly governed by the 
nonlinear geireration of poloidal diamagnetic velocity. In 
a similar approach was taken, without an averaging 
over the sheet, and a value of the velocity at the center 
of the sheath was used. 

To summarize, we have shown that the dynamics of 
magnetic islands can be strongly affected by the presence 
of a background of interchange modes. In the low resis- 



tivity and/or small A' limit, the coupling between the 
magnetic flux and the pressure is dominant compared to 
that between the magnetic flux and the plasma poten- 
tial. In the asymptotic nonlinear regime, a pressure is- 
land structure builds up and the pressure pattern is not a 
flux function except in the center of the magnetic island. 
In fact this regime is a result of a novel nonlinear transi- 
tion that is observed for a wide range of parameters when 
the condition rj < rjc is satisfied. It is initiated by elec- 
trostatic interchange modes which compress the magnetic 
structure and generate small scales inside the island. An- 
other noteworthy finding is that the bifurcation leads to a 
change of symmetry of the diamagnetic velocity that oc- 
curs when the energy of the large scale interchange mode 
is of the same order of magnitude as the thermal energy 
contained in the cell maintaining the magnetic structure. 
The destabilization leads to a poloidal rotation of the is- 
land that is linked to the nonlinearly generated diamag- 
netic velocity in the current sheet. The basic phenomena 
highlighted by our results are reproducible over a large 
region of parametric space and in that sense appear to 
be generic albeit within the constraints of our minimal- 
ist model. Effects ignored in our model including parallel 
heat conduction, parallel ion dynamics and contributions 
of the Hall and electron inertia terms may bring about 
some modifications. Investigation of such effects in an en- 
larged model are therefore necessary to provide a more 
global perspective of this complex phenomenon and for 
which our present studies provide a minimalist and basic 
description. 
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